New Optimization Methods for Converging Perturbative Series with a Field Cutoff 



O 
o 

(N 

Or 
<D 

o 
> 

(N 
(N 
O 

o> 
o 
m 
o 



CD . 
43 ■ 



X 



B. Kessler, L. Li and Y. Meurice * 
Department of Physics and Astronomy, The University of Iowa, Iowa City, Iowa 5224-2, USA 

We take advantage of the fact that in A</> 4 problems, a large field cutoff 4> max makes perturbative 
series converge toward values exponentially close to the exact values, to make optimal choices of 
4>max- For perturbative series terminated at even order, it is in principle possible to adjust (f> max in 
order to obtain the exact result. For perturbative series terminated at odd order, the error can only 
be minimized. It is however possible to introduce a mass shift m 2 — > m 2 (l + r]) in order to obtain the 
exact result. We discuss weak and strong coupling methods to determine <f)max and n. The numerical 
calculations in this article have been performed with a simple integral with one variable. We give 
arguments indicating that the qualitative features observed should extend to quantum mechanics 
and quantum field theory. We found that optimization at even order is more efficient that at odd 
order. We compare our methods with the linear 5-expansion (LDE) (combined with the principle 
of minimal sensitivity) which provides an upper envelope of for the accuracy curves of various Pade 
and Pade-Borel approximants. Our optimization method performs better than the LDE at strong 
and intermediate coupling, but not at weak coupling where it appears less robust and subject to 
further improvements. We also show that it is possible to fix the arbitrary parameter appearing in 
the LDE using the strong coupling expansion, in order to get accuracies comparable to ours. 



I. INTRODUCTION 

Perturbative methods in quantum field theory and 
their graphical representation in terms of Feynman di- 
agrams can be credited for many important physics ac- 
complishments of the 20th century [1,2]. Despite these 
successes, it is also well-known that that perturbative se- 
ries are asymptotic [3,4]. In concrete terms, this means 
that for any fixed coupling, there exists an order K in 
perturbation beyond which higher order terms cease to 
provide a more accurate answer. In practice, this order 
can often be identified by the fact that the K + 1-th 
contribution becomes of the same order or larger than 
the previous ones. The "rule of thumb" consists then in 
dropping all the contribution of order K + 1 and larger, 
allowing errors that are usually slightly smaller than the 
K-ih contribution. 

For low energy processes involving only electromag- 
netic interactions, the rule of thumb would probably be 
satisfactory. On the other hand, when electro-weak or 
strong interactions are turned on, it seems clear that for 
some calculations the errors associated with this proce- 
dure are getting close to the experimental error bars of 
precision test of the standard model [5]. In some cases, 
the situation can be improved by using Pade approxi- 
mants and/or Borel transforms [6]. However such meth- 
ods rarely provide rigorous error bars and do not always 
work well at large coupling or when non-perturbative ef- 
fects are involved. 

In the 21-st century, comparison between precise ex- 
periments and precise calculations may become our only 
window on the physics beyond the standard model. It is 
thus crucial to develop methods that go beyond the rule 



of thumb and provide controllable error bars that can be 
reduced to a level that at least matches the experimental 
error bars. In order to achieve this goal, we need to start 
with examples for which it is possible to obtain accurate 
numerical answers that can be compared with improved 
perturbative methods. This can be achieved with a rea- 
sonable amount of effort in the case of scalar field theory 
(SFT), which we consider as our first target. 

For SFT formulated with the path integral formalism, 
it has been established [7,5] that the large field configu- 
rations are responsible for the asymptotic nature of the 
perturbative series. A simple solution to the problem 
consists in introducing a uniform large field cutoff, in 
other words, restricting the field integral at each site to 
\<S>x\ < 4>max- This yields series converging toward val- 
ues that are exponentially close to the original ones [5] 
provided that (f> max is large enough. Numerical examples 
for three models [5], show that at fixed 4>max, the accu- 
racy of the modified series peaks at some special value 
of the coupling. At fixed coupling, it is possible to find 
an optimal value of (j) max for which the accuracy of the 
modified series is optimal. The determination of this op- 
timal value is the main question discussed in the present 
article. When comparing the three subgraphs of Figs. 2 
and 3 of Ref. [5] which illustrate these features, one is 
struck with the similarity in the patterns observed for 
the three models considered (a simple integral, the an- 
harmonic oscillator and and SFT in 3-dimensions in the 
hierarchical approximation). It is thus reasonable to de- 
velop optimization strategies with the simplest possible 
example, namely the one- variable integral, for which the 
calculation of the coefficients of various expansions does 
not pose serious technical difficulties. As we will see, 
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there exist several ways to proceed and the complicated 
dependence of the accuracy on the coupling constant cer- 
tainly justifies this initial simplification. 

In this article we address the question of the optimal 
choice of the field cutoff <p m ax with the simple integral 
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This integral can be seen as a zero dimensional field the- 
ory. It has been often used to develop new perturbative 
methods [4] , in particular the LDE [8] . The coefficient of 
the quadratic term m 2 is set to I in all the numerical cal- 
culations discussed hereafter, however it will sometimes 
be used as an expansion parameter . 

The effects of a field cut on this integral and the rea- 
son why it makes the perturbative series converge are 
reviewed in section II. Some useful features of the strong 
coupling expansion to be used later are discussed in sec- 
tion III. Our treatment will be different for even and 
odd orders. For series truncated at even orders, the over- 
shooting of the last positive contribution can be used to 
cancel the undershooting effect of the field cut. In other 
words, the errors due to the truncation of the series and 
the field cutoff compensate exactly for a special value 
of the field cutoff ^(A). This value is calculated ap- 
proximately using weak and strong coupling expansion 
in section IV. For series truncated at odd orders (section 
V), the two effects go in the same direction and the error 
can only be minimized. However, an exact cancellation 
can be obtained by using a mass shift m 2 — > m 2 (l + 77). 
We then need to find r)(4>max, A) such that the cancel- 
lation occurs. In practice, it is desirable to have ij as 
small as possible and we will in addition impose that 
dij/d(t) max = 0. This condition fixes the otherwise un- 
specified (pmax- 

The methods presented here have qualitative feature 
that can be compared with the LDE [8] , where the arbi- 
trary parameter can be seen as providing a smooth cut in 
field space, or with variational methods [9] , where weak 
and strong coupling expansions are combined. This is 
discussed in section VI. The main conclusion is that the 
method which consists in determining the value of 4> m ax 
which is optimal for even scries in the weak coupling, 
using the strong coupling expansion provide excellent re- 
sults at moderate and strong coupling. We also show that 
it is possible to fix the arbitrary parameter appearing in 
the LDE using the strong coupling expansion, in order 
to get accuracies comparable to ours. In the conclusions, 
we discuss possible improvement at weak coupling and 
the extension of the model in more general situations. 



The basic question in ordinary perturbation theory is to 
decide for which values of the coupling, the truncated se- 
ries at order K is a good approximation, which in our 
example means 
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The ratios dk+i/dk — — 16k grow linearly when fc — > 00 
and in order to get a good accuracy at order K, we need 
to require A << 1/16K. 

An alternate way of seeing this is that the integrand 
is maximum at <f> = 2\/fc. On the other 
hand, the truncation of e~ Al * at order K is accurate 
provided that A</> 4 << K. The truncated expansion 
of the exponential is a good approximation up to the 
region where the integrand is maximum, provided that 
\{2y/K) 4 «K, which implies A << 1/16K. 

It is useful to represent the above discussion graphi- 
cally. The number of of significant digits as a function of 
the coupling is given in Fig. 1. It is important for the 
reader to get familiar with this kind of graph, because 
we will use them in multiple occasions later in the paper. 
The number of significant digits is minus the log in base 
10 of the relative error. At sufficiently small coupling, the 
behavior becomes linear with a slope which is minus the 
order. Remembering the minus sign above, the intercept 
diminish with the order. It is possible to construct an 
envelope for the curves at various order, in other words, 
a curve that lays above all the curves and is tangent at 
the point of contact. In Fig. 1, we have used a semi- 
empirical formula to draw an approximate envelope: we 
have used the order fc as the (continuous) parameter of 
a parametric curve 



x = -Log 10 (16fc) 

y = -Log 10 [(^(fc+l))-*- 1 



(4) 
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II. EFFECTS OF A FIELD CUTOFF 

In this section, we discuss the effects of a field cut- 
off for the integral defined by Eq. (1). Wc first discuss 
the problems associated with usual perturbation theory. 



A careful examination of the figure at low coupling (e. g., 
near I0~ 2 ) shows that as the order increases, the accu- 
racy increases up to an order where it starts decreasing. 
The envelope is the boundary to a range of accuracy that 
is inaccessible using ordinary perturbation theory. 
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FIG. 1. Number of significant digits versus A at order 1, 2, 
... 15. In all graphs, all the logs are in base 10. As the order 
increases, the curve rotates clockwise and moves left. The 
thick line is the envelope Eq. (4). The large blank area in 
the upper right corner is the region not accessible in regular 
perturbation theory. 



this order, we reach an optimal accuracy. Before doing 
this with a different procedure for even and odd orders 
as explained in the introduction, we will first discuss the 
strong coupling expansion of Eq. (1). 



III. STRONG COUPLING EXPANSION 



In the following, we will often use the strong coupling 
expansion of the integral (1) and a few points should be 
clarified. Our original integral Z(X) vanishes in the limit 
where A — > oo. However X^Z{X) has a finite non-zero 
limit and can be expanded in powers m 2 /A 1 / 2 : 

oo 

A*Z(A)=]T(™ 2 /A 1/2 ))V (7) 



A simple way to convert the asymptotic series into a 
converging one [7,5] consists in restricting the range of 
integration to \cf)\ < cj) ma x- On the restricted domain, 
e~ A< ^ converges uniformly and one can then interchange 
legally the sum and the integral. We are then considering 
a modified problem namely the perturbative evaluation 
of 

/-\-<Pmax 
# e-(™ 2 / 2 ^ 2 -^ 4 . (5) 
-tfimax 

As the order increases, the peak of the integrand of 
a,k (see Eq. (3)) moves across 4> max and the large or- 
der coefficients are suppressed by an inverse factorial: 
| o.fc | < V^(t)m ax /k\. At the same time, we have an ex- 
ponential control of the error: 

/>oo 

\Z(X) - Z(X, <p max )\ < 2e- A ^ / d<f> e-(™ 2 / 2 )^ 2 

(6) 

Everything works in a very similar way for other nu- 
merically solvable A</> 4 problems in D = 1 (anharmonic 
oscillator) and D = 3 (scalar field theory in the hierar- 
chical approximation). The only difference being that in 
these two cases, a more demanding computational effort 
is required. This should be kept in mind while discussing 
the general strategy to be followed. If we could calculate 
as many perturbative coefficients as needed, an obvious 
strategy would be to pick a field cut (j) ma x large enough to 
satisfy some accuracy requirement. Then, given that the 
modified series is convergent, we could calculate enough 
coefficients to get an answer with the required accuracy. 
Unfortunately for any other problem than the integral, 
it is difficult to calculate the coefficients. A more re- 
alistic approach is to assume that we can only reach a 
fixed order and pick the field cutoff in such a way that at 



with 

^ = (-l) i (l/2)'+ 1 (l/Z!)r(Z/2+ 1/4) (8) 

This expansion is converging over the entire complex 
plane. However, if we look at the first few orders dis- 
played later in Fig. 6, one might be tempted to conclude 
that the series has a finite radius of convergence because 
the curves representing the significant digits seem to have 
a "focus" near A = 10 -15 . To be completely specific, 
we mean that in Fig. 6, the four curves labeled SO to 
S4 seem to intersect at a given point. However, as many 
more orders are displayed, the apparent focus moves left 
and a "caustic" (envelope) appears. This is shown in 
Fig. 2. The only difference with Fig. 1 is that the region 
which is inaccessible is now below. In other words, it 
is impossible to reach arbitrarily low accuracy using the 
strong coupling expansion! 



ST. COUPL. INT. 




iog 10 ^ 

FIG. 2. Number of significant digits versus A at order 1, 
11, 21, ... 141 in the strong coupling expansion. As the order 
increases, the curve rotates counterclockwise and moves left. 
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In lattice field theory, the strong coupling expansion 
is similar to the high temperature expansion and for 
D > 2 we expect that the expansion has a finite radius 
of convergence due to the existence of a low temperature 
phase. This difference is not fundamental for the discus- 
sion which follows because we never use the large order 
contributions of the strong coupling expansion. Conse- 
quently, for any practical purpose, the situation will be 
similar to the case where we have a finite radius of con- 
vergence. 



IV. EVEN ORDERS 
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If 4>max is the only adjustable parameter, the perfect 
choice is a solution of: 



FIG. 3. Solutions C { °\D K and E K denned by Eqs. (11), 
(12) and, later in the text, by Eq. (18). 
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k\ 



k r+4> max 



c - (m 2 /2)0 2 (0 4 )fc (1Q) 



Below, we prove that this equation has no solution when 
K is odd, and K is assumed to be even in this section. 
This equation can be solved numerically with good accu- 
racy using Newton's method or a binary search. Our goal 
is to find approximate methods (which can be used in 
more complicated situations) to solve this equation and 
compare them with the accurate numerical solutions. In 
the rest of this section, we consider the cases of strong and 
weak coupling estimates of the optimal value of 4> m ax{ty- 



A. Strong Coupling Estimates 

Multiplying both sides of Eq. (9) by A 1 / 4 and expand- 
ing in powers of m 2 /A 1 ^ 2 we obtain at zeroth order that 
Xfitnax — @fp > with Cjp a solution of 



K 



4£(-i)' 



k=0 



(C (0) )fc +l/4 

fc!(4fc + l) 



r(i/4) . 



(ii) 



The solutions of this transcendental equation are dis- 
played in Fig. 3 for various orders K. Asymptotically, 

~ 0.75 + 0.28K. These solutions can be compared 
with the solutions of the equation 



(D K ) K /K\ , 



(12) 



which can be used as a rough estimate of Xcj)^^ . Asymp- 
totically, D K ^ AK + where A = 0.278465 ... is a 
solution the transcendental equation e~ A ~ 1 = A. 



This lowest order (in the strong coupling) estimate of 
the optimal value of 4>max is quite good. In Fig. 4, we 
see that for K = 6 it provides a significant improvement 
compared to the regular pcrturbative series at order 6 
for A > 10~ 2 . In Fig. 4, we also compare with the ac- 
curacy at fixed cuts. For a fixed value of 4>max, Eq. (9) 
has one solution for a given A and the accuracy becomes 
infinite at this value. In the Fig. we see only peaks of fi- 
nite height, but we see that the approximation goes quite 
high in the peak, in other words, we localize the optimal 
value quite well. 

INT . ORDER 6 
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FIG. 4. Significant digits obtained with the optimal cut 
4>max (A) estimated using a strong coupling expansion at low- 
est order (W6S0), compared to results at three fixed cuts and 
regular perturbation theory (PT6) at order 6. 

We can now proceed to higher orders in m 2 /A 1 / 2 using 
the expansion 



A</> 



max 



EcgV/A 1 ' 9 )', 



(13) 



and plug it in the expansion in the same parameter of Eq. 
(9). The new coefficients obey linear equations which 
can be solved order by order. The optimal (j>rnax 

(A) cal- 
culated at the four lowest orders in m 2 /A 1 / 2 are shown 
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in Fig. 5. As explained in section III, below a certain 
value of A a few orders in the strong coupling expansion 
won't help and one needs much higher order to improve 
the estimate in this region. After a short reflection, one 
can conclude that the "focus" observed in Fig. 6 is com- 
patible with Fig. 5. 

A few words should be said about the notations we use 
for the curves in the figures. When we write W6S1, this 
means that we use the weak coupling expansion up to 
order K=6 in Eq. (9) (this is the W6 part) and a strong 
coupling expansion at order 1 in m 2 /A 1 / 2 (this is the SI 
part) in the calculation of the optimal (j) m ax ■ In addition, 
PT8 means the 8th order in regular perturbation theory. 
In some figures, some of the indexes appear directly near 
the corresponding curve. 

The accuracy of the truncated series at <p ma x calculated 
with the higher order corrections in m /A 1 / 2 in Fig 6. 
For comparison, the accuracy obtained by using only the 
strong coupling expansion Eq. (7). is also shown. The 
figure makes clear that the method proposed here repre- 
sents a significant improvement compared to the separate 
use of the conventional weak and strong coupling expan- 
sions. 



FIG. 6. Significant digits obtained with the optimal cut 
4>max (A) (corresponding to a truncated expansion at order 6 
in the weak coupling) estimated using the strong coupling 
expansion at orders 0, 1, 2 and 3 (solid lines), compared to 
significant digits using only the strong coupling expansion of 
the integral at the same orders in the strong coupling (dashed 
lines) and regular perturbation theory at order 6 (PT6). 



B. Weak coupling 

As we learned in the previous subsection, as A de- 
creases, the optimal value of (f> ma x increases. In the limit 
of a weak coupling, the "tails" of the integral that we 
removed become a small quantity. It is thus advanta- 
geous to split the l.h.s. of Eq. (9) into its bulk and tails 
and expand c~ A< ^ in the bulk where it is justified. The 
resulting (exact) equation for (f> max is then 

/•+oo 00 

2 / dtfe-<™ a / 2 >* 2 -** 4 =- J2 «fc(^-x)A fc , 

J(t> m ax k=K+l 

(14) 
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FIG. 5. Estimates of the optimal </> ma x(A) at the first four 
orders in the strong coupling expansion 
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In the limit of very small A, Eq. (14) becomes 

(2/(mV roax ))e-(™ 2 / 2 ><^ ~ -a K+1 \ K+1 ■ (15) 

The l.h.s. has a functional form similar to semi-classical 
estimates of the energy shifts in quantum mechanics 
[5,10]. A more refined version of this equation is 
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C )A K+1 , (16) 



It is clear that the two above equations have solutions 
only when K is even because in this case cik+i < 0. In 
Appendix A, it is shown that this property extends to 
the exact Eq. (9). More precisely the r.h.s. of Eq. (9) is 
positive for K even and negative for K odd. 

Eq. (16) can be further improved by including higher 
order truncations at odd orders: 



dcj) e 



-(m 2 /2)0 2 -A0 4 



K+3 



ak(<Pmax)^ k , 

(17) 



k=K+l 



and so on. In the following, we refer to the successive 
approximations defined by Eqs. (15), (16) and (17) as 
approximations Al, A2 and A3 respectively. The esti- 
mates of the optimal 4> ma x obtained with these approx- 
imations and the corresponding accuracies as a function 
of the coupling are shown in Figs. 7 and 8. 
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FIG. 7. Estimates of the optimal (f)max(X) obtained with 
the approximations Al, A2 and A3 defined in the text for K 
= 6, compared to numerical values (empty circles). 
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FIG. 8. Significant digits obtained with the approximations 
A\, A2 and A3 for K — 6, compared to order 6 to 9 in regu- 
lar perturbation theory (solid lines rotating clockwise as the 
order increases). 

One can see that Al provides good estimates of <j> m ax 
optimal only at very small A. On the other hand, A2 and 
A3 both provide good estimates even at large coupling. 
Not surprisingly, the accuracy of A2 (A3) merges with 
order 7 (9) in regular perturbation theory. 



V. ODD ORDERS 



As explained in section IV, for series truncated at odd 
K, the r.h.s. of Eq. (9) is negative. The best that we can 
do is to minimize the error (i. e., the difference between 
the r.h.s. and the l.h.s. . The minimization condition 
implies that A^ a:E = Ek, the unique (see Appendix A, 
where Jk is defined) solution of 



m JK{Ek) 



k=0 



k\ 



(18) 



In the following, wc refer to this condition as the Principle 
of Minimal Sensitivity (PMS) condition. This terminol- 
ogy has been used [8] in the LDE where the variational 
parameter is fixed by requiring that the final estimate 
depends as least as possible on this parameter. The solu- 
tions Ek are displayed in Fig. 3 which shows that they 
are asymptotically close to the solutions C R r> _ 1 obtained 
at the lowest order in the strong coupling expansion. The 
accuracy obtained using the PMS condition is by con- 
struction the envelope of the accuracy obtained for all 
possible (f) m ax- This is illustrated in Fig. 9. 
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FIG. 9. Significant digits obtained with the PMS condition 
compared to the same quantity at fixed cuts. 

Nevertheless, an exact match between the original in- 
tegral and the truncated perturbative expansion with a 
field cut can be obtained by using a mass shift m 2 — > 
m 2 (l + rf). Using obvious notations, we denote the 
cut integral defined in Eq. (5) with this mass shift 
Z(X,$ max ,r\). The level curves of the perturbative ex- 
pansion of Z(X, (frmax, v) a t fixed A follow different pat- 
terns at odd and even orders as illustrated in Fig. 10. 

At even order, all the level curves cross the ij = line 
and there is no need for a mass shift. This case was 
discussed in section IV. At odd order, the level curve 
corresponding to the exact value Z(X) defines a curve 
vi&max)- In appendix B, we show that this curve stays 
in the half-plane 77 < 0. We arc now free to pick an ar- 
bitrary value of 4> m ax and adjust r\ — r){<f> max ). In the 
following we will pick <p m ax in such a way that 77 is as 
small as possible. This can be accomplished by solving 
the equation 




dr)/d<j> max = , 



(19) 



for (fimax- In appendix B, we show that this condition is 
indeed equivalent to the PMS condition (18) and conse- 
quently we have simply <j) m ax = 

(Ek/X) 1 / 4 . With this 
choice, the introduction of r] is a natural continuation 
of the optimization at 77 =0. Estimations of 77 for this 
choice of <j) m ax can be obtained approximately at strong 
and weak coupling. 
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FIG. 10. Level curves of the perturbative expansion of 
Z(A,4> max ,r)) at order 6 and 7 in A evaluated at A = 0.1. 



In the limit of arbitrarily small coupling, we can treat 
as a quantity of order X K+1 and use it to make up for 
the "missing" even contribution that would allow a solu- 
tion of Eq. (9). This reasoning implies the weak coupling 
estimate 



V ~ -A K+1 flx + i(^a I )m(2/^) 1/2 



(20) 



On the other hand, at strong coupling, rj grows like A 1 / 2 
and we need to expand 



i 



B" -io- 



-20 



1 r 


1 1 1 1 ■ 


■ 2___ 




0"" 






' 1 




o NUM 




W7S 




W7W8 



-4 -2 2 4 

iog 10 x 

FIG. 11. Optimal i) as a function of A. Numerical val- 
ues (empty circles) compared to approximate values obtained 
with the weak coupling expansion expansion (W7W8) corre- 
sponding to Eq. (20) with K = 7 and the orders 0, 1 and 2 
in the strong coupling expansion (W7S) from Eq. (21). 
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FIG. 12. Significant digits for the same approximations as 
in Fig. 11 compared to the regular perturbative series at or- 
ders 7 and 8 (thin solid line), the PMS result with t?=0 at 
order 7 (small dots) and the envelope of regular perturbation 
theory (thick solid line). 



roV~ 1/2 = I>i(ro 2 M 1/2 )' ■ (21) 



The two approximation work well in their respective 
range of validity as shown in Fig. 11. The significant 
digits obtained with the various procedures are displayed 
in Fig. 12. One sees that the mass shift provides a sig- 
nificant improvement compared to the PMS condition at 
r?=0. If we compare the two methods in their respec- 
tive region of validity, the improvement provided by the 
strong coupling method is more substantial. Not surpris- 
ingly, W7W8 merges with PT8 at weak coupling. 



We can now compare the accuracy of various estimates 
based on strong coupling expansions at the same order in 
m 2 /A 1 / 2 . Examples are shown in Fig. 13 where the accu- 
racy obtained with three methods relying on estimates at 
order one in to 2 /A 1/7 2 are displayed. One can see that as 
the coupling becomes large, the accuracy increases at the 
same rate in the three cases. As we already know in the 
even case, our method significantly improves the basic 
strong coupling expansion from Eq. (7). However, the 
improvement based on even order K = 2q in A performs 
significantly better than the improvement based on the 
odd order K = 2q + 1. Consequently, when in the next 
section we compare with other existing methods, we will 
restrict ourselves to the even case. 
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FIG. 13. Significant digits obtained with the strong cou- 
pling expansion at order one in m 2 /A 1,/2 and the approxima- 
tions W6S1 and W7S1 discussed previously. 
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FIG. 15. Comparison of the delta expansion at order 7 (57) 
with various methods discussed previously and regular pertur- 
bation at order 7 (PT7). 



VI. COMPARISON WITH OTHER METHODS 



There exists several methods to improve the accuracy 
of asymptotic series. These include Padc's approximants 
[6] applied to the series itself or its Borel-transform and 
the LDE [8] . These methods are compared among them- 
selves in Fig. 14. 




FIG. 14. Comparison of the delta expansion at order 7 (57) 
with Pade approximants and the Pade-Borel method. In both 
cases the approximants [4/3], [3,4] and [2,5] have been used. 

One can see that at weak coupling, the LDE provides 
an upper envelope for the accuracy while at strong cou- 
pling it prevails more significantly. Consequently, we only 
need to compare our results to the LDE. This is done 
in Fig. 15 where we see that at strong and moderate 
coupling, our methods provide a significant improvement 
compared to the LDE. On the other hand, at weak cou- 
pling, the improvements that we proposed, do not per- 
form as well as the LDE. 



The results of Fig. 15 have been obtained by making 
the replacements [8] m 2 — > f2 2 + 5(m 2 — O 2 ) and A — > SX. 
We then expanded the perturbative series at order K in 
A to order K in 5. The arbitrary parameter fi 2 was de- 
termined by requiring that the derivative of the estimate 
with respect to fi 2 vanishes. This was called the PMS 
condition in Rcf. [8] and it has a solution at odd orders 
only. 



At even orders, it is however possible to proceed in a 
way similar to what we have done in subsection IV A, 
namely matching the strong coupling expansion of the 
estimate with the usual strong coupling of the integral 
in order to determine the arbitrary parameter. At order 
zero in m 2 /A 1 / 2 , this results into a transcendental equa- 
tion which has a solution at even orders only. Higher 
orders corrections to X/fl can then be calculated by solv- 
ing linear equations just as in subsection IV A. The nu- 
merical results for K — 6 are displayed in Fig. 16. One 
can see that this method and the method presented in 
subsection IV A have very similar accuracy at moderate 
and strong coupling. 

INTEGRAL 
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FIG. 16. Significant digits obtained using the strong cou- 
pling expansion at order one in m 2 /A 1 ^ 2 to determine ft in 
the 5 expansion at order 6 compared the approximation W6S1 
discussed previously. 

The procedure we have used above is closely related to 
variational methods [9] where weak and strong coupling 
expansions were combined for various purposes. The only 
difference is that here we simply imposed the matching 
with the strong coupling expansion rather than resorting 
to cxtremization procedures or large order scaling argu- 
ments. 



VII. CONCLUSIONS 

In conclusion, for even series in the weak coupling, the 
method which consists in determining the optimal value 
of max using the strong coupling expansion provides ex- 
cellent results at moderate at strong coupling. There is 
room for improvement at weak coupling. In particular, 
progress could be made by finding accurate approxima- 
tions to calculate a large number of terms in the r.h.s. of 
Eq. (14). 

The methods used here can be extended to quantum 
mechanics and in particular for the anharmonic oscillator 
where similar calculations have been partially performed 
[14]. We are planning to apply the methods developed 
here for higher dimensional SFT where the LDE seems to 
converge very slowly [11] (see also [12,13] for methods to 
improve the situation). One difficulty is to calculate the 
perturbative coefficients with a field cutoff. Monte Carlo 
methods have been recently developed for this purpose 
[15]. 
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Y. M. was at the Aspen Center for Physics while this 
work was in progress. Y. M. was partly supported by 
a Faculty Scholar Award at The University of Iowa and 
a residential appointment at the Obermann Center for 
Advanced Studies at the University of Iowa, while the 
manuscript was written. 



APPENDIX A: NON-EXISTENCE OF 
SOLUTIONS IN THE ODD CASE 

In this appendix, we show that Eq. (9) has no solu- 
tion when K is odd. For this purpose we introduce the 
truncated exponential series: 

K 

^w^f-ijy/fc!, (ai) 

fc=0 

and their complement 

oo 

g K ( X )= ^(-l) k X k /k\. (A2) 
k=K 



Using the fact that f' K = —fx-i and a similar relation 
for the g, one can show by induction that for K even, fx 
is strictly positive and gx is positive with its only zero at 
zero. For K odd and x > 0, gm is negative and decreases. 
Given that the r.h.s. of Eq. (14) is the integral with a 
positive measure of —gx+i over positive argument, we 
see that the r.h.s. is positive when K is even and nega- 
tive when K is odd. Since the l.h.s. is always positive, 
they are no solutions for K odd. 



APPENDIX B: SPECIAL FEATURES OF r){<f> M Ax) 

The function r){4> max ) is the solution of the equation 
Z(X) = /" e-( m2 /2)(i-H ! (« ra „ ))^ 2 / K (A0 4 ) . (Bl) 

In this appendix, K is assumed to be odd. For x > 0, we 
have Jk{x) < e~ x because gx+i{x) > for K + 1 even 
and x > (see appendix A). We can compensate this un- 
derestimation by making the integration measure more 
positive, in other words picking the parameter rj < 0. 
The l.h.s. of Eq. (Bl) is independent of <j> max . Tak- 
ing the derivative of Eq. (Bl) with respect to (j) ma x and 
imposing that 4> ma x is a solution of dr\jd$ max — , we 
obtain that for this special value of 4> m ax, we have 

fK(X<f> 4 max )e- {m2/2)il+r1 ^ - . (B2) 
which implies the PMS condition (18). 
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